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ABSTRACT 

The growth of galaxies is a key problem in understanding the structure and evolution of the universe. Galaxies 
grow their stellar mass by a combination of star formation and mergers, with a relative importance that is 
redshift dependent. Theoretical models predict quantitatively different contributions from the two channels; 
measuring these from the data is a crucial constraint. Exploiting the Ultra VISTA catalog and a unique sample 
of progenitors of local ultra massive galaxies selected with an abundance matching approach, we quantify the 
role of the two mechanisms from z = 2 to 0. We also compare our results to two independent incarnations 
of semi-analytic models. At all redshifts, progenitors are found in a variety of environments, ranging from 
being isolated to having 5-10 companions with mass ratio at least 1:10 within a projected radius of 500 kpc. 

In models, progenitors have a systematically larger number of companions, entailing a larger mass growth 
for mergers than in observations, at all redshifts. Generally, in both observations and models, the inferred 
and the expected mass growth roughly agree, within the uncertainties. Overall, our analysis confirms the model 
predictions, showing how the growth history of massive galaxies is dominated by in situ star formation at z ^ 2, 
both star-formation and mergers at 1 < z < 2, and by mergers alone at z < 1. Nonetheless, detailed comparisons 
still point out to tensions between the expected mass growth and our results, which might be due to either an 
incorrect progenitors-descendants selection, uncertainties on star formation rate and mass estimates, or the 
adopted assumptions on merger rates. 

Subject headings: cosmology: observations - cosmology: theory - galaxies: general - galaxies: formation - 
galaxies: evolution - galaxies: high redshifts - galaxies: distances and redshifts 


1. INTRODUCTION 

Even though in the last decades much attention has been 
dedicated to the study of galaxy formation and evolution, un¬ 
derstanding when and how the most massive galaxies formed 
and how they evolve with time are still controversial ques¬ 
tions. In the standard paradigm of structure formation, dark 
matter halos assemble hierarchically in a gravitational col¬ 
lapse, and galaxies form inside these structures following the 
radiative cooling of baryons. Stars in today’s most massive 
galaxies ^ IO^^Mq) are formed very early (50% at z ^5, 
80% at z ^3) and in many small galaxies. Model massive 
galaxies can have a number of effective progenitors as high as 
^5 and assemble surprisin gly late. Predictions are model de¬ 
pendent; e . g. acc ording to De Lucia et al. | ( [2006] ); |De Lucia| 
|& Blaizot| ( 2007| ) half their final mass is typically locked-up 
in a single galaxy after z ^0.5. 

Many physical processes have to be taken into account to 
explain the growth of massive galaxies. Star formation is ex¬ 
pected to play an important role at higher redshifts as a large 
fraction o f massive galaxies at z ^ 2 have high star formation 
rates (e.g, [van Dokkum et al ||2004[ |Papovich et al ||2006| ). 
However, the old stellar ages of the most massive early-type 
galaxies (e.g , [Thomas et al |[2005| [van Dokkum & van der 


Marel|2007 ) and the existence of apparently “ red and dead” 


] |2008 [van Dokkum||2008| ) suggest that at least some of the 
growth i s due t o other mechanisms, like m e rgers. Recently, 
Grahamj (|2Q13|); [Pullo & Grahamj ( |2Q13j ); [Graham, Dulloj 


& Savorgnan ( 2015[ ) have also suggested that some massive 
galaxies have evolved by accreting a large disc of gas which 
rapidly forms stars rather than growth only via mergers. 
Below z ^ 1, most massive galaxie s ^ IQ^^Mq) are 


generally found in dense environments ([Blanton & Moustakas 
I, such as at the center of clusters and are identified as the 


galaxies with small sizes at z = 1.5 — 2.5 (e.g, [Cimatti et al 


Brightest Cluster Galaxies (BCGs). Nonetheless, the recent 
survey MASSIVE ( [Ma et a/.~[[2014[ ) showed that this might 
not be always the case. They observed 116 galaxies with 
^ IO^^-^Mq and distance D < 108Mpc (z < 0.025) find¬ 
ing that 48-56% of them are located in groups (35-38% are 
actually central galaxies), while 6-14% are isolated. Many 
studies focused on characterizing the assembly of massive 
galaxi es, both from a theoret i cal and observation al point of 
view. [Lidman et al \ ( [2012[ ); [Lin et al7 ( [2013[ ) found that 
the stellar mass of BCGs increases by a factor of ^ 2 since 
z^ 1. Star formation rates in BCGs at z^l are generally too 
low to result in significant amounts of mass. Instead, most of 
the mass build up occurs through mergers. In semi-analytic 
models, a ccretion of satellite galaxies are mainly dry and mi¬ 
nor (e.g., [De Lucia & Blaizot[[2007[ ), while in observations 
many works point to major mergers in the center of clusters 
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as the main responsible for the mass growth (Rasmussen et al. 


[20101 Brough et al [201 ![ Bildfell et al [2012[ 

Lidman et al 

201 3[). Characterizing three BCGs at z^O.l wit] 
panions, [Brough et al [([2011[) found that the c 
two of the BCGs would merge with the BCG wi 
Additional examples of likely major mergers a 

1 nearby com- 
ompanions of 
Lthin 0.35 Gyr. 
m be found in 

Rasmussen et al [([2010[);[Yamada et al [([2002[); [Collins et al\ 


2009). While it is clear that mergers do occur, it is not yet 


clear what fraction of the stars in the merging galaxies ends 
up in the BCG and what fraction is distributed throughout the 
cluster. High resolution simulations suggest that between 50 
to 80% of the mass of mergers is not locked into galaxies, 
but is distributed throughout the cluste r ( [Conroy, Wechsler ^ 
|Kravtsov|2007|[Puchwein et al. |2010| ). Recent observational 
studies, on the basis of color gradients, support simulations, 
claiming that at least half of the mass is lost into the intraclus- 
ter medium of the clusters (e.g.,|Lidman et al |2012[ [Burke, j 
[Hilton & Collins|2015[ ). 

At higher redshift, the environment in which massive galax¬ 
ies reside is less characterized, and a clear correspondence 
between massive galaxies and BCGs is lacking. Nonethe¬ 
less, a nu mber of studies characteriz ed the build up of massive 
galaxies. [Ownsworth et al. [ (2014), using a variety of num¬ 
ber density selections, claimed that more than half of the total 
stellar mass in massive galaxies {M^ ^ iq11.24 ^^) ^ — q 3 


is created ext ernally to their z = 3 progenitors, [van Dokkum 
\et al \ ( [2Q10[ ), selecting galaxies at a constant number den¬ 
sity of ^ = 2 X \Q~^Mpc~^ , found that the role of mergers 
might be even more important, with star formation account¬ 
ing only for 20% of the total mass growth. Connecting high 
and lo w redshift BCG data v ia evolution of their host halo 
masses, [Shankar et al [ ( [2015] ) found an increase since z ^ 1 
of a factor ^ 2-3 in their mean stellar mass and ^ 2.5-4 in 
their mean effective radius. 

To really understand how individual galaxies have evolved 
and assembled their mass and what mechanisms drive these 
changes, it is important to properly connect today’s most mas¬ 
sive galaxies to their progenitors at earlier times. This re¬ 
quires the non-trivial task of linking galaxies and their de¬ 
scendants/progenitors through cosmic time, which in turn re¬ 
quires assumptions for how galaxies evolve. In recent years, 
a few approaches have been developed to link galaxies across 


sler[[2009| [Behroozi et al 

|2013 

[Wake et al \ 2006 [ [Leja, 

van Dokkum & Eranx 201 3[ 

Munc 

y, Conselice & Ownsworth 


being debated (e.g., Torrey et al. [[2015] ), it is widely rec 


ognized that these approaches are far superior than selecting 
g alaxies at fixed stellar mass for studies of galaxy evolution. 
[Marchesini et al ] ([2014[) assembled the first sample of 


galaxies defined to be the progenitors of galaxies with 
logM^/MQ > 11 . 8 atz = 0 (ultra-massive galaxies, hereafter 
UMGs) from z = 3 using a semi-empirical ap proach based on 
abund ance matc hing in the ACDM paradigm ( [Behroozi et al 


|][TO 1 see & Characterizing the stellar population prop¬ 
erties of the progenitors (masses, ages, dust star formation), 
they claimed that at least half of the stellar content of local 
UMGs was assembled at z > 1, whereas the remaining was 
assembled via merging from z ^ 1 to the present. They also 
found that most of the quenching of the star-forming progen¬ 
itors happened between z = 2.75 and z = 1.25, in good agree¬ 
ment with the typical formation redshift and scatter in age of 
z = 0 UMGs as derived from their fossil records. The progen¬ 
itors of local UMGs, including the star-forming ones, never 


lived on the blue cloud since z = 3. 

Using the unique progen itor-descendant sample presented 
Marchesini et al. [ ( 2Q14[ ), in this paper we focus on the en- 


m 


vironment in which these progenitors reside, and test whether 
their mass growth can be ascribable mainly to mergers or to 
star formation. In particular, we explicitly test the model pre¬ 
dictions for the different contributions to the stellar mass as¬ 
sembly since z ^2. Bo th semi-analytic models (e.g., [Zehavi, 


Patiri & Zheng[[MT^ and abundance matching techniques 


based on halo occupation models (e.g., [Conroy & Wechsler 


2009[ ) indicate that star formation is important at all halo 


masses at z^ 2 , at z < 1 accretion through mergers dominates 
at the high-mass end (^ 10 ^^/z“^M q) of the halo mass distri¬ 
bution, where star formation is negligible, while at intermedi¬ 
ate redshift both contributions are important. Solving possible 
discrepancies found between observations and simulations is 
beyond the scope of this paper, and is deferred to a forthcom¬ 
ing analysis. 

We parameterize the environment in terms of projected dis¬ 
tance from the progenitor, since we have no information about 
the mass of the haloes these galaxies reside. We only consider 
mergers between galaxies with a mass ratio at most of 1 : 10 . 
We also compare our observational resu lts to the predictions 
of two semi-analytic models, name ly the De Lucia & Blaizot| 
( [2007 1 hereafter DLB07) and the [Henriques et al. [(2015| 
hereafter HI5) model, to investigate whether the Marchesini 


\etaL [ ( [ 20 T^ approach to link galaxies across the cosmic time 
is supported by these models. 

Throughout the paper, we assume Hp = 70kms~^M pc~\ 
Uo = 0.3, and Ua = 0.7. We adopt a [Kroupa[ ( [200 1[ ) initial 
mass function (IMF) in the mass range 0.1-100 Mq. 

2. DATASET 

Our sample is drawn from the A^-selected catalog of the 
COSMOS/UltraVISTA field from [Muzzin gF^ (2013a). The 
catalog covers 1.62 deg^ and includes point-spread function- 
matched photometry in 30 photometric bands over the wave¬ 
length range 0.15 /im-24/im from the available GALEX ( [Mar-[ 
[tin et al7~ 


S , Canada-France- Hawaii Telescope/Subaru 
^ _ 07]), Ultra VISTA ([McCracken et f 2 /.~[[ 2012 [ ), 

and S-COSMOS ( [Sanders et al |2007[ ) data sets. Source s are 

selected from the DR l Ultra VISTA X^-band imaging ( [Mc¬ 
Cracken et al [2012 ) which reaches a depth of Ks^tot < 23.4 


at 90% completeness. Details on the photometric catalog con¬ 
struction, photometric redshift measurements, and stellar pop- 
ulation properties’ estimates can be found in [Muzzin et al \ 
(2013 a). Briefly, stellar population properties were derived by 
fitting the observed spectral energy distributions (SEDs) from 
the GALEX UV to the S pitzer-IRAC 8 /im photometry with 
Bruzual & Chariot ( [2003[ ) models assuming exponentially de¬ 


clining SFHs of the form SFR oc , where SFR is the Star 
Formation Rate, t is the time since the onset of star formation, 
and r s ets the timescale of the decline in the S FR, solar metal- 
licity, a [Calzetti et al. [([2000[) dust law , and a [Kroupa[ ( [2001 ) 
I MF (see also [Marchesini et al. |2014[ ). 

[Marchesini et al. ( 2014[ ) investigated the effects of different 
SED-modeling assumptions by adopting, among others, dif¬ 
ferent SFHs, and metallicities. Adopting a delayed-r SFH in 
place of the exponentially-declining SFH allows for increas¬ 
ing SFR at earlier times. The delayed-r model implies SFRs 
smaller by ^ 0.1 dex and stellar ages ^)sfr larger by r ^ 0.1 
dex compared to the default SED- modeling assumptions. Re¬ 
laxing the assumption on metallicity by leaving it as a free 
parameter in the SED modeling does not noticeably change 
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the result, indicating that the impact of fixing the metallicity 
to the solar one is almost negligible. Overall, the systematic 
effects on the stellar population properties are found to be sig¬ 
nificantly smaller than the corresponding typical random un¬ 
certainties for most of the different SED-modeling assump¬ 
tions. Therefore, results are robust and not very sensitive to 
reasonable choices of the SED-modeling assumptions. 

The redshift-de pendent stellar m ass completeness limit has 
been presented in |Muzzin et al \ (2013b). This was deter¬ 
mined by selecting galaxies belonging to the available deeper 
samples and then scaled fluxes and to match the K-band 
completeness limit of the Ultra VISTA sample {Ks^tot = 23.4). 
The upper envelope of points in the {M^^scaied — z) space repre¬ 
sents the most massive galaxies at Kg = 23.4, and so provides 
a redshift-dependent 100% completeness limit for the Ul- 
traVISTA sample. Similarly, [Muzzin et al. \ (2013b) also de¬ 
rived 95% mass-completeness limits for the sample, which 
increases the sample by a factor of 1.4. Given th is substantial 
increase in statistics, we follow [Muzzin et oT] (2013b) and 
adopt the 95% mass-completeness limits. 

The quiescent/star- forming separation has been done us¬ 
ing the rest-fra me U-V versus V-J color -color diagram and 
is pres ented in [Muzzin et al. \ (2013 a), Marchesini et al. \ 


( |2014| ). This method has the ability to separate red galax 
ies that are quiescent from reddened (i.e., dust-obscured) star- 


forming galaxies (see, e.g.. 

Labbe et al [[20061 

Wuyts et al 

20071 [Williams et al |2009t 

Brammer et al |201 It [Patel et al 

1201 lllWhitaker et al. |201 1 

[Muzzin et al [2013a). 


2 . 1 . The selection of the progenitors of local ultra-massive 

galaxies 

The selection of the progenitors of local UMGs has been 
performed adopting a semi-empirical approach which uses 
abundance m atching in the ACDM paradigm (see [Marchesini 
\et al |[2QT^ . This method accounts for mergers and scat¬ 
ter in mass accretion histories. Deta ils on this technique can 
be found in jBehroozi et al. | ( |2013| and references therein). 
Briefly, the galaxy cumulative number density at redshift z\ is 
converted to a halo mass with equal cumulative number den¬ 
sity using peak halo mass functions. Then, for halos at that 
mass at zi, the masses of the most-massive progenitor halos 
at Z2 > Zi are recorded using to the halos’ mass accretion his¬ 
tories. Einally, the median halo progenitor mass at zi is con¬ 
verted back into cumulative number densities using the halo 
mass function at Z 2 - 

The progenitors of the low-z population of very massive 
galaxies are traced by identifying, at each redshift, the stel¬ 
lar mass for which the evolving cumulative number density 
intersects the cumulative number densi ty curves deriv ed from 
the Ultra VISTA stellar mass functions ( [Muzzin et al. [ 2013 a). 
In this way, a sample of progenitors of galaxies with a mass 
of - at z ^ 0 is assembled. 

The typical error on the progenitors’ stellar mass resulting 
from the uncertainties of the observed stellar mass functions 
and cumulativ e number densities was fo und to be in the range 
0.03-0.07 dex ( [Marchesini et al |2014[ ). The inferred growth 
in stellar mass of the progenitors was therefore found to be 
0.45±0.13 dex and 0.27±0.08 dex from z = 2 and z=l, re¬ 
spectively, to z = 0. If the scatter in mass accretion histories 
is also included in the error analysis, the uncertainties on the 
inferred growth in stellar mass of the progenitors increase by 
a factor of —1.7 ( [Marchesini et al |2014[ ). 

2 . 2 . Our sample 



Fig. 1.— Redshift histogram of our progenitors of UMGs sample (black). 
The distribution of quiescent (red) and star-forming (blue) galaxies is also 
shown. 


In this work, we use the sample of progenitors of UMGs 
defined in Marchesini et al. | (2014), but we limit our analy¬ 
sis to galaxies at z <2. In this way, at all redshifts, the Ul¬ 
tra VISTA sample includes all galaxies at least as massive as 
1:10 the mass of the closest progenitor in the redshift range 
±0.05 X (1 ^ Zpr ) with Zpr the photo-z of the progenitor. The 
redshift range is chosen to take into account the typ ical photo- 
z accuracy , which is redshift dependent (see, e.g., [Muzzin g7[ 
al [[2013K 191 progenitors enter our selection. As shown 


in Eig. [IJ the number of progenitors depends on redshift: 
there are 11 galaxies at 0.2 < z < 0.5, 41 at 0.5 < z < 1, 
69 at 1 < z < 1.5 and 70 at 1.5 < z < 2. The different fre¬ 
quency is mainly due to the different volume covered by the 
different redshift bins, which is — 10 x larger at z ^ 1.75 
than at z ^ 0.35. Indeed, the volume probed is —25 Gpc^ 
at 0.2 < z < 0.5, —125 Gpc^ at 0.5 < z < 1, —200 Gpc^ at 
1 < z < 1.5, and —235 Gpc^ at 1.5 < z < 2. 

As already pointed out by [Marchesini et al \ ( [2014| ), for 
z < 1 all galaxies are quiescent, while at higher redshift the 
fraction of star-forming progenitors is not negligible. 

To characterize the environment in which progenitors are 
embedded, we also make use of the entire Ultra VISTA cata¬ 
log, to which we apply a mass cut to ensure 95% of complete¬ 
ness (see 

Hereafter, we will refer to those galaxies around the progen¬ 
itors within a specified projected radius, at least as massive 
as 1:10 the mass of the progenitor and in the redshift range 
±0.05 X (1 -\- Zpr ) as companion galaxies. 

We note that there might be systematic effects in the data 
that alter the robustness of the results, such as systematic er¬ 
rors in photometric redshifts and contamination of the pho¬ 
tometry from emission lines. The latter might result in over¬ 
estimates of stellar masses. [Marchesini et al [ ( [2014[ ) already 
investigated the possible systematics in the Ultra VISTA sam¬ 
ple, showing that at the redshifts he re considered th ey should 
not impact on our findings (see also [Muzzin et al [ 2013a). 


3. THEORETICAL PREDICTIONS 

In this paper we compare the findings from the observations 
with the predictions from theoretical models on the environ¬ 
ment of the progenitors of today’s UMGs and on the relative 
importance of merging vs. in-situ star formation to their in¬ 
ferred growth in stellar mass to investigate the physical pro¬ 
cesses implemented in the models. To this aim, the model 
predictions are derived in two different ways. In one case, 
the model-predicted assembly histories of the progenitors are 
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directly exploited to determine the overall growth in stellar 
mass of the descendants. In the other case, the growth in stel¬ 
lar mass is obtained following the same assumptions used in 
the observations. 

We exploit galaxy catalogs from t wo semi-analytic mod - 
els run on the Millennium Simulation ( [Springel et al. [2005] ). 
This uses 10^^ particles of mass 8.6 x lO^Tr^^liT^Tolracettie 
evolution of the matter distribution in a cubic region of the 
Universe of 500/z“^M/7c on a side from z = 127 until z = 0, 
and has a spatial resolution of 5 h~^ kpc. 

We use two different semi-analytic models to investigate 
how different assumptions about the physical processes acting 
on the baryonic component impact the evolution of the galaxy 
masses. 

The semi-analytic model discussed in DLB 07 builds on the 
metho d ology and prescriptions introduced in 
( |2001| ); [be Lucia et al | ( [2QQ4| ); [Croton et al. 


Springel et al. 
( |200^ and has 


been the first variant of the “Munich” models family that has 
been made publicly avai lable. The DLB07 m odel is based 
on WMAPl cosmology ( [Sanchez et al |[2006] ) and includes 
prescriptions for supernova-driven winds, follows the growth 
of supermassive black holes and includes a phenomenologi¬ 
cal description of AGN feedback. The model neglects envi¬ 
ronmental physical processes such as ram pressure and ha¬ 
rassment, but assumes that when galaxies are accreted on to 
a more massive system, the associated hot gas reservoir is 
stripped instantaneously. This induces a very rapid decline 
of the star formation histories of satellite galaxies, and con¬ 
tributes to create an excess of re d and passive galaxie s with 
respect to the observations (e.g., [Wang et al [2007] ). The 


DLB07 model is mainly tuned to reproduce the K-band lu 
minosity function at z=0. 

We also use the model presented in HI5, which represents 
one of the most recent updates of the Munich models. The 
HI5 model uses the Planck first-year cosmology and basically 
contains the same physics as the DLB07 model but has a more 
sophisticated treatment for the evolution of satellites. Differ¬ 
ently from the DLB07 model, which does not include a chan¬ 
nel for ICL formation, the HI5 model includes tidal stripping 
as a channel for ICL. In addition, it adds a modification of 
the time-scale to re-accrete gas ejected through galactic winds 
and modifies the ram-pressure stripping in halos less massive 
than ^ IO^'^Mq. The model is tuned to reproduce recent data 
on the abundance and passive fractions of galaxies and the 
galaxy stellar mass function from z = 3 down to z = 0. We 
refer to the origina l papers for more details . _ 

As e xplained in Springel et al ( 2001| ); De Lucia et al 
( |2004| ), models make a distinction between centrals, satellites 
and orphans. Centrals (type 0) are located at the position of 
the most bound particle in their halo. These galaxies are fed 
by gas cooling from the surrounding hot halo medium. Satel¬ 
lites (type 1) were previously central galaxies of another halo, 
which then merged to form the larger system in which they 
currently reside. For these galaxies gas is no longer able to 
cool on to halo galaxies. Orphans (type 2) are galaxies no 
longer associated with distinct dark matter substructures, and 
in the DLB07 model their stellar mass is not affected by the 
tidal stripping that reduces the mass of their parent halos. In 
the HI5 model such orphans are unable to retain gas ejected 
by supernova feedback, which is moved to the hot halo of the 
galaxy group. Tidal forces can completely disrupt the stellar 
and cold gas components of orphan galaxies, which are then 
added to the intra-cluster light and the hot gas atmosphere of 
the group/cluster, respectively. In both models, orphans may 



Fig. 2.— Predicted median progenitor mass growth by semi analytic mod¬ 
els (left: DLB07, right: H15) for galaxies in the same mass range of our 
progenitors. Masses at z = 0 are the masses of the descendants. Error bars 
on the y-axis represent the 25^^ and 75^^ percentiles. At each redshift, 10 
random extractions have been performed. 

later merge into the central galaxy of their halo. In our analy¬ 
sis, when useful, we will distinguish among the three types of 
galaxies. 

3.1. Our sample 

For both the DLB07 and the H15 models, we extract from 
the available catalogs all the galaxies at z = 0.36, 0.76, 1.28 
and 1.77 with stellar mass in the same mass range spanned by 
the progenitors at the corresponding redshift (approximately 
within ±0.15 the median mass of the progenitors). 

Since observed masses might be characterized by system¬ 
atic errors, from the DLB07 model we also extracted samples 
of galaxies to test the impact of these errors. We assigned 
to each galaxy mass in the models a random Gaussian error 
with width 0.03 x (1 ±z) (following jllbert a/. |2013| and 
then considered only those galaxies whose perturbed mass is 
in the mass range ±0.15 the stellar mass of the progenitors at 
the corresponding redshift. We performed the entire analysis 
using both samples, without finding noticeable differences be¬ 
tween the results. Therefore, in the following, we will present 
only the analysis performed on the sample with the original 
masses from the models. 

Overall, in the DLB07 (H15) model 1027 (1076) galaxies 
have been extracted at z = 0.36, 608 (1231) at z = 0.76, 447 
(1969) at z= 1.28 and 311 (1732) at z = 1.77. 

From the models, we also extract the information regarding 
the merger trees and the descendants of these galaxies down 
to z ^ 0, with the aim to investigate the real mass growth pre¬ 
dicted by the models. In addition, we also get the virial mass 
of the haloes in which these galaxies reside and those of their 
descendants, to further characterize from a theoretical point 
of view the progenitors’ hosting environment. 

Finally, we select all galaxies within a box of 1 physical 
Mpc on a side, centered on each massive galaxy considered, in 
order to characterize also the environment of the progenitors 
of UMGs in the same way as in observations. 

4. DESCENDANTS OF THE PROGENITORS IN THE MODELS 

First of all, we can tes t whether the approach adopted by 
Marchesini et al ( |2014| ) to link galaxies across the cosmic 
time is supported by the models. 

From simulations, we randomly extract the same number of 
progenitors found in observations at the corresponding red¬ 
shift and compute the median mass of their descendants at 
z ^ 0. We repeat the sampling ten times, to take into account 
sample variance. Figure shows that the selection based on 
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TABLE 1 

Percentages oe type 0, type 1 and type 2 galaxies among the progenitors and descendants in the two semi-analytic models. 


DLB07 


) progenitors 


descendants 


H15 


progenitors 


) descendants 


0.36 

0.76 

1.28 

1.77 


type 0 

type 1 

type 2 

type 0 

type 1 

type 2 

type 0 

type 1 

type 2 

type 0 

type 1 

type 2 

9ll[ 
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12 li 

9 0 + 0.8 

00 

1+1 
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1 + 


0 5^^'^ 

87^2 
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87l[ 

lolj 

9 q+0.8 
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10 l[ 

4 3+^'^ 

93^2 


0 9+^-^ 
^•^-0.5 

85 I 3 

y +2 

-2 

8+2 

88 l[ 

0 1 -|-0.9 
^•^- 0.8 


84l[ 

12 l[ 

4 5+^-^ 
^•^- 0.6 

95tl 


0 6 “^^'^ 
^•^-0.5 

851' 

9+3 

^-2 


90t[ 

6 

'^•^- 0.8 

0 1 -|- 0.6 
■^•^- 0.6 

00 

1+1 

4 6 +^-^ 


the abundance matching method does indeed select galaxies 
whose mass evolution is consistent with what expected from 
the DLB07 model: progenitors at the different redshifts will 
turn into galaxies at 0 with masses of r\j IqII-^Mq, that is 
the mass inferred by [Marchesini et aL~ ( 2014 ). The model 
also reproduces the inferred stellar mass at the intermediate 
redshifts. In contrast, in the HI5 model, the median mass of 
the descendants at z = 0 is lower than ^ 10 ^ ^-^Mq in most of 
the extractions, even though within a large dispersion. It is in¬ 
teresting to note that the abundance matching method does not 
seem to work for the HI5 model, which instead is the model 
that should be in better agreement with the statistics of the 
halo occupation model by construction (see H15 for details). 

Differences between the two models might be due to the 
fact that the H15 has introduced tidal stripping and therefore 
has a smaller number of satellites, producing a smaller mass 
growth through mergers. 

In the following sections we will quantify the separate role 
of star formation and mergers in the galaxy mass growth, also 
from a purely theoretical point of view, by explicitly inspect¬ 
ing the merger trees of a subsample of galaxies. 

Table shows the percentage of galaxies of a given type, 
for all objects extracted from the simulations. In the DLB07 
model, at all redshifts, the vast majority (>90%) of the pro¬ 
genitors are type 0, with the fraction slightly decreasing going 
from higher to lower redshifts. This finding is probably due 
to the fact that at lower redshift the number of satellites is 
larger and that these galaxies had more time to grow. It sug¬ 
gests that while at higher redshift massive galaxies are most 
likely at the center of their halo, in the local universe there is 
a larger fraction of massive galaxies that are satellites. Only 
< 1% are type 2. Among the descendants, the fraction of type 
0 is smaller and decreases with increasing redshift of the pro¬ 
genitors, ranging from 85% to 88%. In contrast, the fraction 
of type 2 galaxies is much larger and they are as common as 
satellites (^ 7% at all redshifts). This suggests that the while 
almost all massive galaxies that are central at z = 0.36 will be 
central also at z = 0, the probability of central galaxies to turn 
into satellites increases with increasing redshift. 

The HI5 model presents a systematically lower fraction of 
type 0 galaxies (85-90%) and a slightly larger fraction of type 
1 among the progenitors, and a similar fraction of type 0 and 
higher of type 1 among the descendants. 

Figure shows the halo mass distribution of all the se¬ 
lected progenitors and also the mass distribution of the haloes 
where the descendants of the progenitors reside, for both mod¬ 
els. Both in the DLB07 and HI5 prescriptions, the typical 
mass of the haloes hosting the progenitors clearly increases 
with decreasing redshifts: galaxies at z = 1.77 are found in 
haloes of mass ^ 10^^-^ — for the DLB07 model, 

and of ^ 10^^-^ — 10 ^‘^-^Mq for the H15 model; galaxies at 
z = 0.36 are found in haloes of mass ^ — IO^^-^Mq 


DLB07 


DLB07 


12.5 13.0 13.5 14.0 14.5 15.0 15.5 

log(JlWA7^) 


12.5 13.0 13.5 14.0 14.5 15.0 15.5 
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Fig. 3.— Normalized halo mass distribution for progenitors at different 
redshfits (left panels) and for their descendants at z = 0 (right panels), as 
indicated in the labels. Upper panels: DLB07 model, bottom panels: HI5 
model. 


for the DLB07 model and of ^ 10^^ — 10^^ for the H15 
model; supporting the finding that the growth of the ma ssive 
galaxies is coup led to the growth of their haloes (e.g, |Tin- 
|ker et al |2012| and references therein). Looking at the halo 
mass distributions for the descendants, in the DLB07 model, 
they span a similar halo mass range, even though the peaks of 
the distributions slightly depend on the redshift of the pro¬ 
genitors: descendants of the z=1.77 progenitors are found 
in slightly less massive haloes than the descendants of the 
z=0.36 progenitors (medians values are M^ir = 
and Mv/r = respectively). Therefore, even though 

the progenitors turn into massive galaxies of similar mass, 
they actually might not end up in the very same galaxies, as 
expected by the fact that the growth of structure is stochas¬ 
tic. It also might suggest that also the environment around 
galaxies should be taken into account when linking galaxies 
across time, and not only the stellar mass. In contrast, in the 
HI5 model all the descendants span a similar halo mass range 
(Mv/r ^ 10^^ — 10^^-^) with a median halo mass of 
Even though in the HI5 model not all the progenitors will end 
up in massive galaxies as selected by the abundance match¬ 
ing technique, these distribution show that they will end up in 
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very similar environments. 

5. THE ENVIRONMENT AROUND THE PROGENITORS OE UMGs 

In order to understand the processes that induce the ob¬ 
served galaxy mass growth, we first characterize the envi¬ 
ronment in which the progenitors of UMGs reside. If these 
galaxies are found in overdense regions, they should easily 
undergo mergers; on the other hand, if they are most likely 
isolated, their growth should be attributable to other factors, 
such as in-situ star formation. 

Figure]^ shows some examples of false color images in the 
BzK filters covering a FoV of 500 kpc on a side of galaxies re¬ 
siding in different environments at different redshifts. In order 
to demonstrate the range of environments of these galaxies, 
for each redshift bin, we selected a galaxy with no other com¬ 
panions within a projected sphere of 250 kpc, a galaxy with 


In Figure are also overplotted the same quantities for data 
drawn from the models. The upper panel shows the results 
for the DLB07 model, the lower panel for the H15 one. At 
each redshift, we randomly extract the same number of pro¬ 
genitors found in observations at the corresponding redshift, 
and then we compute the projected distances^ to define the 
number of companions within a certain projected radius. We 
repeat the sampling ten times, to take into account sample 
variance. We then consider the maximum range spanned by 
the extractions. Qualitatively, within a large spread, both the 
DLB07 and H15 models follow the observational trends, at all 
distances. Nonetheless, in both models there are extractions 
where the number of galaxies is systematically higher than 
observed, at all distances, especially in the DLB07 one. This 
might be related to the well known issue of the over-prediction 


3-4 companions a galaxy with ^8-10 companions. Clearly, 

mann et al 

2011J and has been found to be present also in 

progenitors reside in a variety of environments, which will 

clusters (Vu 

cani et al [[2014|). Discrepancies are solved by 

have a different role in their growth throughout the cosmic 

construction in the HI5 model. 


time. 

5.1. The total number of satellites around progenitors 

Using the observed sample drawn from the Ultra VISTA 
catalog, we compute the number of companions around each 
progenitor at different redshifts. We consider portions of sky 
centered on the progenitor and of different physical radii and 
count the number of companions, namely galaxies with mass 
ratio >1:10 and redshift within the range ±0.05 x (1 ^Zpr), 
that fall into the projected area. 

Figure shows the results considering galaxies within 50, 
100, 250 a nd 500 kpc from the progenit or, respectively. Ac¬ 
cording to [Munoz-Cuartas et al. ( [2011| ), haloes with M^ir ^ 
IO^^Mq have typically a virial radius of ^ 0.5 Mpc below 
z ^ 1, and of ^ 0.35 Mpc at z = 2, while haloes with M^ir ^ 
IO^^Mq have typically a virial radius <0.15 Mpc at all red¬ 
shifts. Our binning has been chosen to inspect regions of sky 
that correspond to the virialized region around the progenitors 
at the different redshift, assuming they are located in a vari¬ 
ety of haloes. The smallest distance indeed should sample the 
virial radius for any of our progenitor galaxy, while the largest 
should sample well beyond it. 

To get rid of the different volumes covered at the differ¬ 
ent redshifts and therefore the different number of progeni¬ 
tors, for each sample values are normalized to the total num¬ 
ber of progenitors at the considered redshift. In observations 
(thick lines in both panels), at any redshift, ~80% of progen¬ 
itors have no galaxies closer than 50 kpc. At z < 0.5 there 
might be an excess of progenitors with 3 galaxies within 50 
kpc 10%). All these progenitors have been found to live 
in X-ray selected COSMOS poups with halo m asses in the 
range 10^^ — lO^'^M 2 ooc /^0 ( George et al [2011 ). Enlarging 
the radius of interest, the number of companions around pro¬ 
genitors increases and a dependence on redshift might appear 
for the most extreme galaxies with the largest number of com¬ 
panions. Going from higher to lower redshift, distributions 
shift toward a larger number of companions, suggesting that 
the environment around progenitors gets richer. Considering 
our largest radius (500 kpc), we find that about 25% of pro¬ 
genitors at 1.5<z<2 have no companions and another 25% 
have at most two companions. On the other hand, at lower 
redshift at most 5% of progenitors are isolated. At z < 0.9, 
^45% of progenitors with at least one companion have been 
found to live in X -ray selected COSMOS groups ( [George ^Tj 
al [2011[ ). No group catalogs are available at higher redshift. 


The fact that in simulations there is a non-negligible spread 
among the different extractions indicates that the sample vari¬ 
ance is not marginal, therefore a larger sample of observed 
galaxies will be needed to draw more robus t conclusions. 

Our results are in line with the results of [Tal et ari ( [2013] ), 
who found that the total number of galaxies within a mass 
range of 1:10 and within roughly 400 kpc of the massive 


galaxy is on average 2-3 in all redshift bins (see al so [Tal et 
\al. |2Q12[ Quills & Trujillo|2012[[Gobat et al [2015 1. 


5.2. Comparisons with the stellar mass function 

To further assess the environment in which the observed 
progenitors live, we estimate how many satellites progenitors 
should have based on the stellar mass function down to 1:10 
the mass of the progenitor. This number could give us hints 
whether most of them live in clusters/groups or in the field. 
TST et al \ ( [2014[ ) presented the stellar mass function of 


satellites around central galaxies of different stellar masses 
at 0.2 < z < 1.2. They identified central galaxy candidates 
from the Ultra VISTA catalog. Galaxies were considered to 
be central if no other, more massive, galaxies could be found 
within two projected virial radii. Virial radius estimates at 
a given stellar mass a nd redshift were d etermined using the 
semi-analytic model of [Guo et alTl^lOl l| . They found that the 
mass distribution of satellite galaxies is independent of red¬ 
shift for any given value of central galaxy mass. If we make 
the hypothesis that our progen itors are the centra l galaxies of 
a group, we can integrate the Tal et a/.~[ ( 2Q14| ) mass func¬ 
tion from the mass of the progenitor down to 1:10 of its mass 
and compare the expected number to our observed ones. An¬ 
alytically, we find that we should expect 2-3 companions per 
central galaxies within two p rojected virial radii.^ Following 
[Munoz-Cuartas et ar] ( [201 1[ ), at z ^ 1, the typical virial radius 
of halos with Myir ^ IO^^Mq is ^0.5 Mpc. If we therefore 
consider a radius of 1 Mpc, we find that for z > 1 15% of the 
progenitors have at most one companion, while for z < 1 there 
are no isolated progenitors. This might suggest that most of 
progenitors live in massive structures like groups or clusters, 
but at least at high redshift there is a non-negligible fraction 
of them which have no companions. 

At z < 1 we can push down the mass limit to 1:100, without 

^ For each galaxies we consider projected distances both on the xy, xz and 

planes. __ 

^ Two virial radii is the size chosen by Tal et al. j2014| . 
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Fig. 4.— Examples of BzK images for isolated galaxies within 250 kpc (first column), galaxies with few companions (second column) and galaxies with a 
large number of companions, in four bins of redshifts. Companions are galaxies around the progenitors, at least as massive as ~1:10 the mass of the progenitor 
and in the redshift range ±0.05 x (1 -\-Zpr)- The Ultra VISTA id, the redhisft and the number of companions is indicated in the labels. The FoV of each thumbnail 
corresponds to a radius of 250 kpc. 


being biased by mass-incompleteness. With this mass thresh¬ 
old, analytically, we should expect ~ 9 companions per cen¬ 
tral galaxy within two projected virial radii. In our observa¬ 
tions we find an average of 50±5 galaxies within the same 
radius. This might suggest that at least our low redshift pro¬ 
genitors are located within structures largely dominat ed by 
the presence of small galaxies. As we will see later on ( |6.3| ), 
even though there is not much mass enclosed in these galax¬ 
ies, they play a role in the mass growth of the progenitors, 
given their high SSFR. 

5.3. Number of satellites as a function of distance 

We now focus our attention on progenitors with at least one 
companion within 500 kpc and investigate, on average, the 
variation of the number of companions per projected volume 


with distance, both in observations and in models (Fig. |^. In 
observations, trends with redshift are not detected. However, 
if we sum up the number of companions, we find 4.6±0.9, 
4.8±0.5, 3.7±0.3, 3.6±0.3 galaxies within 500 kpc per pro¬ 
genitor from 0.2<z<0.5 to 1.5<z<2, respectively. Trends 
with distance are also detected. The number of companions 
drops between distances of 100 and 200 kpc and is about con¬ 
stant at larger distances. 

Models show stronger trends with distance. In the DLB07 
model (upper panel of Fig|^ there is on average a larger 
number of companions around each progenitor than in obser¬ 
vations, especially at small distances (<100 kpc), where the 
number of galaxies is more than a factor of two larger than 
in observations. Considering 10 different random extractions. 
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Fig. 5.— Normalized distribution of the number of companions (=galaxies at least as massive as 1:10 the mass of the progenitor and with z within ±0.05 x 
(1 ±Zpr)) at different redshifts and within different radii (as indicated on the top of each panel). Thick dashed lines represent observations; thin solid lines and 
shaded areas represent models (upper panels: DLB07, bottom panels: H15). Numbers in parenthesis give the number of progenitors in each redshift bin. 


the median number of galaxies around progenitors within 500 
kpc is 6±1, 7.4±0.5, 7.3±0.4, 5.9±0.3 from z = 0.36 to 
z = 1.77, respectively. This again might be due to the over¬ 
prediction of the number of satellites in the DLB07 prescrip¬ 
tion. In contrast, in HI5 (bottom panel of Fig |^, the number 
of objects is typically consistent with the observed one (me¬ 
dian numbers are 3.5±0.8, 3.7±0.4, 3.6±0.3, 3.1±0.2). In 
both models no trends with redshift are detected, even though 
the DLB07 model shows a possible inversion in the lowest 
redshift bin. 

Our findings are consistent with the idea that progenitors 
are indeed centrals (so surrounded by a satellite population 
with some decreasing number density profile). Indeed, as 
showed in Table [T] in both models and at all redshifts, >90% 
of the selected galaxies are classified as centrals in their halos. 

The discrepancies found in the DLB07 model are no¬ 
ticeably alleviated if we exclude type 2 galaxies (plots not 
shown), which are mostly found at close distances from the 
center of the halos. Similarly, also in the HI5 model the num¬ 
ber of galaxies at ^50 kpc is strongly reduced when type 2 
galaxies are removed. We note, however, that type 2 galax¬ 
ies cannot be excluded: it has been shown that their presence 
is fundamental to reproduce well several properties, i.e. the 
clustering at small scales of the structures and the differences 
between the galaxy and subhalo profiles in the inner regions 


L different 


of clusters (e.g. Gao et al |2004t|Wang et al. |2006| ). 

We can also investigate whether progenitors with dif 
star-forming properties live in different environments. Mod¬ 
els do not provide us with the necessary information to dis¬ 
tinguish between star-forming as quiescent galaxies as it was 
done for observations, nonetheless, both in models and ob¬ 
servations we can use galaxy Specific Star Formation Rate 
(SSFR = SFR/M^) to distinguish between star-forming and 
quiescent galaxies. We assume that galaxies with \ogSSFR < 
— 11 are quiescent, being this the minimum of the distribu¬ 
tion. In observations, this number roughly corresponds to the 
adopted U-V vs. V-J cut. 

In observations, as already mentioned, while at z < 1 all 
progenitors are quiescent, at higher redshift 60% of progen¬ 
itors are star-forming. Also in models the fraction of quies¬ 
cent galaxies depends on redshift, spanning from ^95(78)% 
to ^ 30(17)% in the DLB07 (HI5) model going from low 
to high-z. Both in observations and simulations, there is no 
significant difference in the trends shown in Fig. [^when we 
select only quiescent candidate progenitors. 

Similarly, we now focus on the properties of the compan¬ 
ions, checking whether their distribution around progenitors 
depends on their star-forming properties. Figure [7] shows the 
variation with distance and redshift of the number of quies¬ 
cent galaxies around quiescent progenitors both in observa- 
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Fig. 6. — Number of companions per progenitor per kpc^ as a function 
of distance at different redshifts, as indicated in the labels. Only progenitors 
with at least one companion within 500 kpc, galaxies at least as massive as 
1:10 the mass of the progenitor and with z within ±0.05 x (1 Zpr) are 
considered. Errors are poissonian. A horizontal shift is applied to the points 
for the sake of clarity. Thick dashed lines and points represent observations; 
thin solid lines and shaded areas represent models (upper panel: DLB07, 
bottom panel: H15). In models, error bars represent the range spanned by the 
10 extractions. 

tions and simulations. Trends resemble those found for the 
total population. In observations the distance dependence is 
less steep, while in models there are hints it might be steeper. 
Both models show an over abundance of satellites at small 
distances. Being overall the trends similar in Fig. |^an d0 we 
can conclude that there are no evident signs of clustering in 
observations, where the quiescent and star forming galaxies 
are similarly distributed around progenitors, while in mod¬ 
els it seems that quiescent galaxies might be more clustered 
around quiescent progenitors. In observations, similar results 
have been obtained when using a U-V vs. V-J cut, showing 
that the results are not much sensitive to the cut adopted to 
separate star forming from quiescent galaxies. 

To summarize, we have found that the number of compan¬ 
ions around progenitors does not depend on redshift. In obser¬ 
vations, going from higher to lower redshift and from smaller 
to larger distances, the environment gets proportionally richer 
of galaxies. Nonetheless, there is a fraction of progenitors that 
do not have companions, suggesting that their mass growth is 
hardly related to merger events. Excluding isolated progen¬ 
itors, we found that the distribution of companions per pro¬ 
jected volume is almost independent on redshift. 

In models, the fraction of isolated progenitors is much 
lower, indicating progenitors live in denser environments. 

6 . THE DRIVERS OE THE PROGENITORS’ MASS GROWTH 


Eig. 7.— Number of quiescent companions per quiescent progenitor per 
kpc^ as a function of distance at different redshifts, as indicated in the labels. 
Points, lines, error bars and colors are as in Fig.[^ 


In this section we aim to investigate which are the most im¬ 
portant factors that drive the galaxy mass growth for progeni¬ 
tors from z ^ 2 to z ^ 0. We will first focus only on mergers 
( |6.1| ), then we will quantify the contribution of in -situ star 
formation ( |6.2| ) and finally combine the two ( |6.3| ), to esti¬ 
mate their relative importance at the different redshifts. 


6.1. What fraction of galaxy mass growth is due to mergers? 

Kitzbichler & White] ( [2Q08| ) investigated the major merger 
rates using catalogs based on the DLB07 model to obtain the 
characteristic time-scale needed by two galaxies of a given 
mass ratio and redshift to merge on the basis of their pro¬ 
jected distance. In the model, to determine whether or not 
two galaxies might merge, it is assumed that when the sub¬ 
halo which hosts a galaxy is tidally disrupted near the center 
of a more massive halo, the galaxy becomes eligible to merge 
with the central galaxy of that halo. Nonetheless, the merger 
does not occur immediately, but rather after a “dynamical fric¬ 
tion time” estimated from the relative orbit of the two objects 
at the moment of subhalo disruption. 

Forz < 1, stellar masses above 5 x IO^Mq and samples lim¬ 
ited to radial velocity difference Av < 3000km/sf 


(Tinerge) = 3.2Gyr- 




50kpc V4.6-1010 Mo 


-0.3 


(‘ + 20 > 


^ The authors suggest to use the timescales for Av < 3000 km/s when 
analyzing data from photometric redshift samples, since the “background” 
correction will not eliminate physically associated galaxies at large velocity 
separation. 
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Fig. 8 .— Observed median progenitor mass growth at different redshifts, as indicated in the labels, assuming that galaxies within 50, 100, 250, 500 kpc will 
merge by z = 0 onto the progenitor. The mass at z = 0 is the sum of the entire mass (stars and dashed regions) or half of the mass (squares and shaded regions) 
of all galaxies with mass ratio 1:10 and redshift ±0.1. Error bars on the x-axis represent the wi dth of the redshift bin, en -or bars on the y-axis represent the 25^^ 
and 75^^ percentiles. The dotted horizontal lines represent the mass of the UMGs at z ~ 0 (from |Marchesini 12014) . 


where Tmerge is the time scale, Vp the projected physical sep¬ 
aration, the stellar mass of the pairs, z the redshift of the 
p rogenitor. _ 

[Kitzbichler & White] ( [2008| ) do not provide a formula 
for higher redshift galaxies, therefore we use the same 
parametrization also at 1< z <2. Given the range of masses 
and redshifts in our sample, we find that galaxies located more 
than ~350 kpc apart at z ^ 2 should not go through a merger 
event by z = 0. 

We can therefore estimate the mass growth due to mergers 
by summing the stellar mass in companion galaxies that are 
expected to merge by z = 0 with the mass of the progenitor."^ 
We analyze two cases: in the first one we assume that the en¬ 
tire mass in the companions will end up in the mass of the pro¬ 
genitor, in the second case only half of it, while the other half 
is assumed to go into some diffuse component. These ar e two 
extreme cases that should bracket the real situation (e.g.,|Con 


roy, Wechsler & Kravtsov|2007t|Lidman et al. |2013[|Burke, 

Hilton & Collins|2015D . 

Figure shows the median mass growth that observed pro¬ 
genitors are expected to have from their redshift to z~ 0.^ 
In the local Uni verse, the typical stelly mass of UMGs is 
^ 10^(Marchesini et al. 2014|. 

We consider different bins of distance, and therefore only 
galaxies within a certain distance from the progenitor. Going 
to very large distances allows us to be as inclusive as possible 
in terms of companions to count and give an estimate of how 
much the mass growth changes as a function of distance. 

First of all, we note that the dispersion (described as the 
25^^ and 75^^ percentile of the distributions) around the me¬ 
dian mass is asymmetric, and is mainly due to the fact that 
we are on the exponential tail of the mass function, therefore 
distributions are not normal. 

The figure shows that the mass growth depends on the con¬ 
sidered distance. Taking into account only galaxies within 
50 kpc from the progenitor, the mass growth is negligible for 
galaxies at z >0.5. In contrast, mergers alone might explain 
the growth of galaxies from 0.2<z<0.5toz = 0. Increas- 


^ For each pair progenitor- companion we estimate the time they need to 
merge and evaluate if they will merge, eventually. 

^ We note that our values slightly differ from those presented in|March-| 
lesini et al |d2014l because they adopt mean masses, whereas we use median 
values. 


ing the radius of interest, the mass growth due to mergers in¬ 
creases; nonetheless it is generally still insufficient to justify 
the expected mass growth. This is true both assuming that the 
entire mass in the companions will end up in the mass of the 
progenitor, and that only half of it will. 

Taking into account all galaxies that can actually merge 
with the progenitor from their redshift to z = 0, on average, 
galaxies increase their mass of 41±9%, 33±1%, 27.9±0.4%, 
31.4±0.6% from z^ 0.35, 0.75, 1.25, 1.75 respectively. 

We note that only for z < 1 we can include in the computa¬ 
tion galaxies with smaller mass ratio (down to 1:100) with re¬ 
spect to the mass of the progenitor, without being affected by 
sample incompleteness. Considering also these galaxies in the 
computation does not strongly influence the results (plots not 
shown), simply because despite there are many more galax¬ 
ies, their low-mass is negligible with respect the mass of the 
progenitor. 

When inspecting models (Fig. [9]), similar results are ob¬ 
tained for the lowest redshift bin and the mass growth for 
galaxies at z = 0.36 is compatible with the mass of the local 
UMGs. In addition, the DLB07 model can explain the mass 
growth in terms of mergers from z = 0.76 (1.28) to z = 0 when 
a radius >100 (250) kpc is considered, the H15 one when a ra¬ 
dius >250 kpc is considered. We emphasize the large spread 
that characterizes the models, which indicates the variety of 
growth histories that characterizes galaxies. 

Recall that in the left panel of Fig. where we compared 
the stellar mass of the progenitors and their descendants as 
given by the models, we found that in the H15 model the mass 
growth is such that the descendants at z = 0 do not have mass 
of ^ IO^^-^Mq, indicating that in this model the progenitor 
selection did not work properly. This is due to the fact that in 
the HI5 model not all the mass in the merging galaxies ends 
up into the central one, but a not negligible fraction goes into 
the intracluster medium (HI5). 

In summary, the analysis above suggests that mergers alone 
can explain the mass growth of the progenitors of the local 
UMGs only for galaxies at z < 0.5 in observations, and z < 1, 
in simulations. In contrast, they do not produce enough mass 
for galaxies at higher redshift. This is due to the fact that at 
z > 1 progenitors are isolated or have very low-mass com¬ 
panions, whose mass is not sufficient to explain the expected 
trends. These results suggest that other factors might play an 
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Fig. 9.— Predicted median progenitor mass growth at different redshifts, as indicated in the labels, assuming that galaxies within 50, 100, 250, 500 kpc will 
merge by z = 0 onto the UMG. The mass at z = 0 is the sum of the entire mass (stars and dashed regions) or half of the mass (squares and shaded regions) of 
all galaxies with mass ratio 1:10. Error bars on the y-axis represent the maximum the 25^^ and 75^^ percentiles of the 10 representations. Upper panel: DLB07, 
bottom panel: HI5. The dotted horizontal lines represent the mass of the UMGs at z ~ 0 (from |Marchesini et al. |2014) . Black triangles represent the inferred 
mass at z = 0 from observations, assuming that all the mass of the satellites fall into the progenitor. 
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important role in the galaxy mass growth. 


6.2. What fraction of galaxy mass growth is due to star 
formation? 


Figure 10 shows the amount of mass progenitors are ex¬ 
pected to gain for star formation, both in observations and 
models. In observations, we adopt both the SFR estimates 
obtained from the SED fitting and those obtained from the 
UV-fIR. The latter are systematically higher than those ob¬ 
tained from the SED fitting and give us an upper limit of the 
growth. In simulations we use the SFR estimates provided by 
the two semi-analytic models. 

We consider a constant SFR in the range of time between 
the redshift of the galaxy and the lowest limit of the next red- 
shift bin. In this way our estimates most likely represent an 
upper limit of the real situation.^ In the computation, we take 
into account the fact that the stellar mass of a galaxy changes 
with time also simply due to the evolution of its stars: as 
they progressively evolve and eventually di e, they retain only 
part of their mass as remnant. Following JPoggianti et al \ 
( |2013| ), who used the |Bruzual & Charlotj ( |2003| ) model, the 


fraction of initial stellar mass that remains is equal to 1 for 
ages less than 1.9 x 10^ yr, while it can be approximated as 
f{t) = 1.749-0.124 X log^ at older ages, where t is the age 
of the stellar population in years. Approximately, in 0.6 Gyr 
galaxies retain ~60-70% of the mass they have formed. 

The left panel of Fig. shows the results for obser¬ 
vations. Estimates from SED fitting produce a little mass 
growth, while estimates from UV-fIR can explain the mass 
growth from one redshift bin to the next. In the DEB 07 
model (central panel of Eig. the median values obtained 
are comparable to those observed from SER estimated from 
the SED fitting. As in observations, star formation alone can 
only marginally explain the mass growth. In addition, results 
of the 10 different extractions are quite similar, as indicated 
by the moderately small scatter. This suggests that in the 
DLB07 model galaxies of similar mass have similar SER at 
the time they have been selected. In contrast, values in the 
HI5 model (right panel) are systematically larger than the ob¬ 
served ones, even though there are some extractions where 
they show agreement. Note that in the HI5 model the scatter 
is very large, indicating that galaxies of similar mass can have 
a wide range of SERs. 


6.3. Combining the contribution of SFR and mergers 

In the previous subsections we have found that neither 
mergers nor star formation alone are able to fully explain the 
expected progenitors’ mass growth from ^ = 2 to z = 0, both in 
observations and, to some extent, in simulations. Here we aim 
to test whether the combined contribution of in-situ and envi¬ 
ronmental processes can produce the expected growth. We 
also consider the mass growth due to star formation in the 
galaxies that will merge with the progenitors. We note that 
our analysis does not consider the contribution from starburst 
during mergers. This is probably not a dominant channel for 
mass growth, but might play a somewhat more important role 
at higher redshift. 

As in the previous section, we compute the mass growth 
due to star formation in the time interval between the redshift 
of the galaxy and the next redshift bin. In this case we use the 

^ In observations, considering a declining SFR does not strongly change 
the results (plots not shown). 


same time interval also to estimate the contribution of merg¬ 
ers, so that we can sum them up together. The left panel of 
EigurepT] shows the results for the observations. Considering 
all galaxies that might eventually merge, when the SER deter¬ 
mined from the SED fitting is adopted, the combination of the 
two contributions marginally explains the mass growth, barely 
tracing the lower limit of the growth. Instead, when we adopt 
the SER estimates obtained from the combination of UV and 
IR luminosities, we recover the expected mass growth. Recall 
that the SFRuv+ir represents an upper limit of the true values, 
given that at these redshifts the AGN contamination might not 
be negligible. However, the SER estimated from the UV-i-IR 
is arguably less biased against heavily obscured star forma¬ 
tion. Therefore, the real growth is expected to be bracketed 
between these two cases. 

The central and right panels of the same Eigure show the re¬ 
sults for the two semi-analytic models. Both predict a growth 
sufficient to support the expected mass growth, at all redshifts. 
Again, the large scatter that characterizes the models suggests 
that the different extractions we performed from the catalogs 
can give quite different results. 

The result obtained in the right panel for the HI5 model is 
in disagreement with the results presented for the same model 
in Eig|^ When we apply the same prescriptions as done in 
the observations to account for the growth from both merging 
and in-situ star formation from one redshift bin to the follow¬ 
ing one, the inferred mass is in agreement with the mass of the 
progenitors at that redshift. As a consequence, we can grow 
all progenitors to a mass of ^ at z = 0. In con¬ 

trast, when we directly consider the mass of the descendants 
at z = 0 as provided by the models, we find a systematically 
smaller mass for the galaxies of the same initial mass. This 
means that some of the assumptions made to estimate these 
contributions may not be sufficient or may be incorrect (e.g., 
stripping and/or merger rates) when adopted for this particular 
model, as it will be discussed in the following section. 


7. DISCUSSION 


The main results of our analysis are summarized in Eigures 


12 and 13 The former presents the relative contribution of 


star formation and mergers to the mass growth, showing the 
inverse of ratio of the progenitor stellar mass at a given red¬ 
shift to the inferred stellar mass at the same redshift. The 
expected stellar mass growth from the abundance matching 
technique, the total measured stellar mass growth obtained 
considering mergers and star formation and the separate con¬ 
tribution of star formation and mergers are shown. The latter 
figure shows the cumulative mass growth due to the two con¬ 
tributions, separately. In both figures, the quoted uncertainties 
represent the maximum and minimum growth, while errors 


on the medians are estimated as 1.253a/^/N, where a is the 
standard deviation ab out the median and N is the number of 
galaxies ( |Rider|1960| ). 

In observations (left panels) the total mass growth we ob¬ 
tained is below the expectations, when the SER values are 
obtained from the SED fitting. Trends are driven by high-z 
galaxies: from 1 <z< 1.5toz^0.5we measure a growth of 
^ 10% while the expected growth is of ^ 40%. This entails 
that at z = 0 ^ 20% of mass is lacking. In contrast, discrep¬ 
ancies are largely reduced when SERs are measured from a 
combination of UV and IR luminosities. In this case, the in¬ 
ferred mass growth is even larger than the expected one at the 
highest redshift. 

Eocusing on models (central and right panels), both pre- 
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Fig. 10.— Median progenitor mass growth due to a constant star formation at different redshifts, as indicated in the labels. Left panel: observations. Stars and 
dashed regions show the mass growth adopting SFR estimates from SED fitting; squares and shaded regions show the mass growth adopting SFR estimates from 
UV+IR. Central panel: DLB07 model, right panel: H15 model. In the panels showing the models, black triangles represent the values obtained from observations. 
Error bars on the x-axis represent the width of the redshift bin (only in observations), error bars on the y-axis represent the 25^^ and 75^^ percentiles. 




Fig. 11.— Median progenitor mass growth due to both star formation and merger at different redshifts, as indicated in the labels. Also the satellites’ mass 
growth is taken into account. A constant SFR and the assumption that all the mass of the satellites will end up in into the mass of the progenitor have been 
made. Left panel: observations. Stars and dashed regions show the mass growth adopting SFR estimates from SED fitting; squares and shaded regions show the 
mass growth adopting SFR estimates from UV+IR, central panel: DLB07 model, right panel: H15 model. In simulations, black triangles represent the values 
obtained fro observations. Error bars on the x-axis represent the width of the redshift bin (only in observations), error bars on the y-axis represent the 25^^ and 
75^^ percentiles. 


scriptions are able to fully explain the mass growth as pre¬ 
dicted by the abundance matching technique, and, possibly, 
even over-predicting it. 

Investigating separately the contribution of star formation 
and mergers, we find that they play a very different role at the 
different redshifts. In observations, the average star formation 
rate is similar to the net growth rate at z = 1.5 — 2 but signif¬ 
icantly smaller at later times. Progenitor’s star formation is 
only important at the highest redshifts, where it might be able 
to explain alone all the mass growth. At lower redshifts merg¬ 
ers acquire importance, and they are the major responsible of 
the observed evolution at z < 0.5. At the intermediate red¬ 
shifts, the growth can be explained advocating the combined 
contribution of star formation and mergers. 


Overal l, these findings are in line with many other stud- 


ies (e.g., van Dokkum et al ||1999 

van 

9okkum||2005l |van 

Dolium et al ||2010t |Tran et al 

12005 

Bell et al \ 

2006; 

White et al |2007l McIntosh et al 

2008 

Naab et al \ 

2007, 

Naab, Johansson & Ostriker|2009[ Ownsworth et al |2014|), 


even though some works have suggested that major mergers 
may play a more prominent role with up to ^ 60% of a mas¬ 
sive galaxies stellar mass growth at z < 2 arising from major 
merger events (e.g., |L6pez-Sa njua n et al ||2012t IFerreras et 


al |2014t|Ruiz, Trujillo & Mdrmol-Queralt6|2014 ). 


In both models the contribution of star formation to the total 
mass growth decreases with time. However, it plays a larger 
role in the HI5 model than in the DLB07, at all redshifts, 
and it shows a steeper decline with time. It goes from 20% 
to 5% in the H15 model, from < 10% to 0% in the DLB07 
model. In the DLB07 model the contribution of mergers in¬ 
creases with time, ranging from 20% to 50%, with a slope that 
is similar to the observed one. However, differently from ob¬ 
servations, they are much more important than star formation 
even at higher redshift. In contrast, in the HI5 model the con¬ 
tribution of mergers is roughly constant with time, showing a 
bump only in the lowest redshift bin. Values are similar to the 
observed ones. 

As already mentioned, the non-negligible spread measured 
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Fig. 12.— Percentage of the mass growth as a function of time as expected by the abundance matching technique (black stars), as measured considering 
the combined contribution of SFR, mergers and SFR in the satellites that will merge (blue points), as measured considering only SFR in the progenitors (green 
squares) and only mergers (red diamonds), for observations (left panel), the DLB07 model (central panel) and the H15 model (right panel). In observations, 
values obtained both considering the SFRs from the SED htting (filled circles) and those from the UV+IR (empty circles). A shift has been app lied to the points 
for the sake of clearness. Error bars represent the maximum and minimum growth, obtained propagating the errors on the medians jRid er p^ . 






Fig. 13.— Mass build- up over time due to star formation and mergers. Upper panels: observations; left panel: SFR from the SED fitting, right panel: SFR 
from the UV+IR luminosities. Bottom panels: simulations; left panel: DLB07, right panel: H15. Black stars and solid lines represent the expected mass growth 
from|Marchesini et al. U2014J with the uncertainties; red regions the contribution of mergers to the mass growth, green regions the contribution of SFR to the 
mass growth. T hick dashed l ines represent the upper limit of the contributions, while thick dash-dotted line the lower limit, both obtained propagating the errors 
on the medians ^Rider|1960) . 
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among the different extractions in simulations indicates that 
the sample variance is not marginal, therefore a larger sam¬ 
ple of observed galaxies will be needed to draw more robust 
conclusions. 

In models, we can explicitly investigate the separate role of 
mergers and star formation by inspecting the merger trees of a 
subsample of galaxies in the highest redshift bin. We consider 
only progenitors at z = 1.77, since most of the progenitors 
at lower redshift are actually descendants of these galaxies, 
therefore will enter their merger trees later on. 

In both models, for each progenitor we select all the galax¬ 
ies in the merger tree that have its same direct descendant in 
the next snapshot and we sum up the masses of the merg¬ 
ing galaxies to compute the mass growth due to mergers. We 
then compare the mass of the merged galaxies to the mass of 
their unique descendant and assume that the difference is due 
to star formation.^ Note that in the HI5 model we consider 
only the galaxies that indeed merged, and not those disrupted 
before merging onto its descendant, whose matter went into 
the intracluster medium. We then start from the descendant 
and repeat the loop down to z = 0, in order to trace the entire 
growth. 

Figure shows the median relative mass growth and the 
median total mass growth for 70 galaxies^ in the DLB07 and 
HI5 models, respectively. As found in the previous sections, 
the contribution of mergers and star formation is different in 
the two prescriptions, being mergers relatively more impor¬ 
tant in the DLB07 model than in the H15 one. In HI5 star 
formation plays an important role also at low redshift. Most 
importantly, we find that in the two models galaxies are char¬ 
acterized by an overall very different mass growth. In line 
with the results shown in Fig. from z=1.77toz = 0 
galaxies in the DLB07 model growth x 1.5 more than in the 
HI5 model. This is likely due to the modifications in the 
HI5 model, which add the tidal stripping and reduce the mass 
of the merging satellites, therefore producing a smaller mass 
growth. Nonetheless, this result is quite surprising given the 
fact that the model reproduces the evolution of the stellar mass 
function and has been calibrated to be in better agreement 
with the halo occupation distribution results (HI5), which 
the abundance matching technique relies on. Understanding 
the reasons of these discrepancies is beyond the scope of this 
work, and is deferred to a forthcoming work. 

7.1. Some caveats 

Overall, our analysis reveals some tensions between the 
mass growth expected by the abundance matching technique 
and the mass growth measured taking into account in-situ star 
formation and mergers, both in simulations and in observa¬ 
tions. The largest discrepancy is seen between 1 < z < 1.5 and 
0.5 < z < 1 in observations, when the SFRs from the SED fit¬ 
ting are used. We note that these two contiguous redshift bins 
bracket a break in the galaxy property distribution: at z > 1 
many progenitors and satellites are still star forming, while at 
z < 1 all the progenitors and most of t he satellites are quies¬ 
cent (see also [Marchesini et al ||2014| ). As we will see later 
on, this transition regime might be responsible for this gap. 

Overall, discrepancies might be due to a number of fac¬ 
tors. First of all, it might be that the abundance matching 
technique adopted to link galaxies across time does not work 

^ Note that this means the stars formed during star-bursts associated with 
mergers are going to be in the mass growth phase 

^ This is the number of the observed progenitors at 1.5< z <2. 


properly. Anyway, our analysis shows that candidates are se¬ 
lected well, at least for DLB07 model. Quite surprisingly, the 
HI5 does not support the selection via the abundance match¬ 
ing technique. Indeed, the mass growth estimated for the HI5 
model and the one that is intrinsic in the model (as shown in 
Figures 13 and 2) are not in agreement. Understanding the 
weakness of the selection criteria is beyond the scope of this 
paper, since the method has been largely discussed in the lit¬ 
erature and it has been found to provide an excellent match 
to a number of galaxy clustering statistics at multiple epochs 


([Kravtsov et al. |2004||Tasitsiomi etal. |2004t|Vale & Qstriker 


pro 


onroy, Wechsler & Kravtsov|2006[JBerrier et al. 

Marin et al. |2008[]Trujillo^Gomez et al. |2Q1 Ij) and to 


a number of population p^e rties (e.g., [Conroy & Wechsler 

2009t[Drory & Alvarez|2008[ ). 


Alternatively, discrepancies might be due to some assump¬ 
tions made. For example, it might be that in reality merg¬ 
ers play a less important role at highe r redshift than that esti¬ 
mated by [Kitzbichler & White[ ( [2008] ) at z < 1. Even though 
merger rates are not expected to vary much with redshift (e.g.. 


Guo & Whitel2008l[Kitzbichler & White|2008l [Wetzel, Cohn[ 


& White[[2009[ ), we might be over-estimating the number of 

mergers. 

In addition, in our treatment, we are not considering some 
other factors that indeed might play a role. The most impor¬ 
tant is the contribution of galaxies whose mass is lower that 
1:10 the mass of the progenitor and therefore do not enter our 
selection. These galaxies can be characterized by high SSFR 
values, hence they double their mass rapidly, therefore giving 
a non-negligible contribution to the total growth. We checked 
that at least at z < 1, where our sample is not affected by in¬ 
completeness, including in the computation all galaxies with 
a mass ratio of 1:100 better reconcile the expected to the ob¬ 
served growth (plot not shown). We can not extend to higher 
redshift because of incompleteness effects. We note, however, 
that at z > 1 the contribution of satellites with mass ratio larger 
than 1:10 might play a more important role than that at lower 
redshift. Indeed, most of them are star forming and are prob¬ 
ably characterized by high SSFR values, therefore giving a 
large contribution to the total growth. At z < 1 many satellites 
are quiescent and contribute less to the total growth. 

Additionally, mergers can also induce bursty events of star 
formation, which can pump up galaxy masses. However, it 
is very hard to properly model these bursts and quantify their 
role in the overall galaxy growth. Not considering the con¬ 
tribution from starburst during mergers has a larger impact 
at z > 1 that at z < 1. Indeed, at higher redshift mergers 
most likely involve star forming galaxies and are accompa¬ 
nied by bursts that enhance galaxy star formation, with a con¬ 
sequent larger mass growth, while at lower redshift, mergers 
most likely take place between quiescent galaxies, therefore 
bursts are rare. 

Finally, in observations, uncertainties in the star formation 
histories, dust content and distribution, the IMF, and other ef¬ 
fects can easily introduce systematic errors of a factor of ^2 
in the star formation rates, particularly at high redshift (see. 


e.g., [Reddy et al \ 

20081 [Wuyts et al [[20091 Muzzin et al 

2009J. Even thoug 

h nowadays there is reasonable agreement 


between the global stellar mass density inferred at any partic¬ 
ular time and the time integral of all the preceding instanta¬ 
neous star-formation activity, modest offsets may still point 
toward systematic uncerta inties that are not negligible (see 
[Madau & Dickinson|2014| for a review). 

Moreover, the tensions seen between 1 < z < 1.5 and 0.5 < 
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Fig. 14.— Predicted mass build-up over time due to star formation and mergers, as obtained from the merger trees, for the DLB (left panels) and H15 (right 
panels) models. Upper panels: cumulative mass fraction, lower panels: cumulative mass growth. The median of 70 merger trees is shown (solid lines), along 
with 1(7 dispersion (dashed lines). Black lines: total growth, red lines and areas: growth due to mergers, green lines and areas: growth due to star formation. See 
text for details. 


z < 1 in observations are at least partly due to the fact that 
our analysis relies on COSMOS data, which covers only one 
field of view. Therefore, we are not able to control for sample 
variance. [Guzzo et al \ ( |2007| ) identified a large-scale struc¬ 
ture at z ^0.73, which certainly contaminates the counts at 
0.5 < z < 1, most probably having an impact on the mass 
functions and on the cumulative number densities involved 
in the selection of the progenitors. Having a larger sample of 
galaxies, based on several fields, is mandatory to really prove 
the existence of the observed gap and to understand its origin. 


8 . SUMMARY AND CONCLUSIONS 

The aim of this paper was to test the model predictions for 
the different contributions to the stellar mass assembly since 
z ^2, and investigate the role of the star formation and merg¬ 
ers at the different redshifts. We compared observational re¬ 
sults with the data of two different semi-analytic models, to 
obtain a better insight on the physical processes responsible 
for the evolution. 

First, we characterized the environment of the progenitors 
of local UMGs at 0.2 < z < 2, selected with a semi-empirical 
approach using abundance matching in the ACD M paradigm 
( [Behroozi et al. ||2Q13t [Marchesini et al. ||2014| ). We inves¬ 
tigated the number of companions around each progenitor, 
in order to give an estimate of the environment surrounding 
these massive galaxies. The number of galaxies with mass at 
least 1:10 the mass of the progenitor and with redshift within 
±0.05 X (1 ^Zpr) around progenitors depends on distance. In 
observations, at any redshift ^ 80% of the progenitors have 
no galaxies within a projected radius of 50 kpc. This number 
drops to 25% at z ^ 1.75 and 5% at lower redshift when a ra¬ 
dius of 500 kpc is considered. In general, going from higher 


to lower redshift the environment gets proportionally richer 
of companions. Models qualitatively agree with observations, 
even though the fraction of isolated progenitors is much lower 
in the models, indicating progenitors live in denser environ¬ 
ments, pointing to the well known over-estimation of satellites 
at high redshift (in the DLB07 model). 

Considering only progenitors with at least one companion 
within 500 kpc, in both observations and simulations the num¬ 
ber of companions decreases with distance. Nonetheless, the 
DLB07 model overestimates the number of companions al- 
most at all distances (in agreement with pr evious results, e.g. 
Weinmann et al. 2011[ Vulcani et al. 2014). The star-forming 


properties of progenitors and companions seem not to influ¬ 
ence the trends in observations, while in models the fraction 
of quiescent companions might be higher around quiescent 
progenitors. 

In the second part of the paper we investigated which are the 
most important factors in the progenitors’ mass growth at the 
different redshifts, characterizing the separate contribution of 
star formation and mergers. 

Overall, our analysis confirms the model predictions, show¬ 
ing how the growth history of massive galaxies is dominated 
by in situ star formation at z ^ 2, both star-formation and 
mergers at 1 < z < 2, and by mergers alone at z < L Nonethe¬ 
less, detailed comparisons reveal some tension between the 
mass growth expected by the abundance matching technique 
and that measured both in observations and in simulations. In 
observations we recover a systematically smaller mass growth 
when SFRs from SED fitting are adopted, whereas we obtain 
an overall comparable mass growth when SFRs from UV-i-IR 
are used. The true mass growth has to be bracketed between 
these two cases. In models, both the prescriptions explain the 
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mass growth as predicted by the abundance matching tech¬ 
nique, and, when errors are taken into account even over¬ 
predict it. The role of the different contributions is different in 
the two prescriptions, highlighting how much the mass growth 
is model dependent. It is worth noticing that the mass growth 
estimated for the HI5 model and the one that is intrinsic in the 
model are not in agreement. This implies that at least some of 
the assumptions made to estimate the different contribution to 
mass growth might be wrong. 

Discrepancies might be due to a number of factors, such as 
an incorrect progenitors-descendants selection, an underesti¬ 
mate of minor mergers (>1:10), the adopted assumptions on 
merger rates, or uncertainties on star formation rate and mass 
estimates. 

In the future, a larger sample of observed galaxies will be 
needed to draw more robust conclusions. Indeed, the non- 
negligible spread measured among the different extractions in 
simulations indicates that the sample variance is not marginal. 
A larger sample would also allow to better investigate the 
impact of the environmental processes on galaxy evolution. 
It has been shown that the efficiency and the time scale of 
quenching of star formation in satellites is halo mass depen¬ 
dent, therefore at any redshift the role of star f ormation is cer- 
tainly different in different environments (e.g.,|Dekel & Bim-| 
|boiin|2006i[Dekel et al 12009] ). 

Another natural step forward to this analysis will be to bet¬ 
ter characterize the gas content and the interstellar medium 
(ISM) in the progenitors’ population, especially in the star¬ 
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forming progenitors, during a period in cosmic history that is 
most critical for the formation of their stars (1 < z < 3). Ob¬ 
servation e.g. with ALMA will allow to carefully investigate 
the ISM, which is the crucial ingredient fueling the activities 
of star formation and AGNs and separating the two contribu¬ 
tion will better constrain the actual role of star formation in 
the total mass growth. 
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